Dissipation effects in random transverse-field Ising chains 
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We study the effects of Ohmic, super-Ohmic, and sub-Ohmic dissipation on the zero-temperature 
quantum phase transition in the random transverse-field Ising chain by means of an (asymptotically 
exact) analytical strong-disorder renormalization-group approach. We find that Ohmic damping 
destabilizes the infinite-randomness critical point and the associated quantum Griffiths singularities 
of the dissipationless system. The quantum dynamics of large magnetic clusters freezes completely 
which destroys the sharp phase transition by smearing. The effects of sub-Ohmic dissipation are 
similar and also lead to a smeared transition. In contrast, super-Ohmic damping is an irrelevant 
perturbation; the critical behavior is thus identical to that of the dissipationless system. We discuss 
the resulting phase diagrams, the behavior of various observables, and the implications to higher 
dimensions and experiments. 
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I. INTRODUCTION 

Continuous phase transitions display the remarkable 
feature of universality: The physics sufficiently close to 
the transition point is independent of microscopic details, 
it only depends on a small number of key parameters such 
as the symmetry of the order parameter, the dimension- 
ality of the system, and the presence or absence of frus- 
tration (for reviews see, e.g., Refs. 1,2). For this reason, 
simple prototypical models are extensively used in the- 
oretical studies of continuous phase transitions as their 
critical behavior can be expected to exactly reproduce 
that of experimental systems. 

Because realistic systems often contain considerable 
amounts of quenched (time-independent) disorder (ran- 
domness), it is important to establish whether or not such 
disorder is relevant or irrelevant at a continuous phase 
transition. In other words, is disorder one of the unim- 
portant microscopic details or is it one of the key param- 
eters that determine the critical behavior? Interestingly, 
the answer to this question is not unique but depends 
on the transition at hand. Early insight was gained from 
analyzing the fate of the disorder strength under coarse 
graining: According to the Harris criterion,^ disorder is 
perturbatively relevant at a critical point if its clean cor- 
relation length exponent ly fulfills the inequality dv < 2 
with d the space dimensionality. In this case, the crit- 
ical behavior of the disordered system must differ from 
that of the clean one. However, the ultimate fate of the 
transition can not be inferred from the Harris criterion. 

In recent years, it has become increasingly clear that 
rare, strongly interacting spatial regions play an im- 
portant role in phase transitions in disordered systems. 
These regions can be (locally) in the ordered phase while 
the bulk system is still in the disordered phase. Their 
dynamics becomes very slow because it involves coher- 
ently changing the order parameter in a large volume. 
Griffiths^'''' showed that these rare regions lead to non- 
analyticities in the free energy not just at the critical 



point but in an entire parameter region around the tran- 
sition which is now known as the Griffiths phase. At 
generic thermal (classical) transitions, the resulting Grif- 
fiths singularities in thermodynamic quantities are very 
weak and probably unobservable in experiment.^ 

Rare regions are more important for zero-temperature 
quantum phase transitions, the critical behavior of which 
is determined by order-parameter fluctuations in space 
and tinre. As quenched disorder is perfectly correlated in 
the time direction, disorder effects are enhanced. The re- 
sulting strong quantum Griffiths singularities^^^ of ther- 
modynamic quantities can take power-law form. The 
critical point itself remains sharp but can be of exotic 
infinite-randomness type as was shown by Fisher. ^'''^^ 

At quantum phase transitions, statics and dynamics 
are intimately coupled. Therefore, changes in the quan- 
tum dynamics can cause changes in the thermodynamic 
behavior. For example, dissipation, i.e., damping of the 
order-parameter fluctuations due to additional degrees of 
freedom, can further enhance the rare-region effects. In 
systems with Ising symmetry, each rare region acts as a 
two-level system. In the presence of Ohmic dissipation, 
it undergoes the localization quantum phase transition of 
the spin-boson problem. ^^ Thus, the quantum dynamics 
of sufficiently large rare regions freezes, ^^'^^ leading to a 
smearing of the global phase transition. ^^ 

In order to gain a more complete understanding of 
the non-perturbative physics of these dissipative rare re- 
gions, the heuristic arguments of Refs. 13-15 need to be 
complemented by an explicit calculation of a prototypi- 
cal microscopic model. Schehr and Rieger made signif- 
icant progress in this direction by applying a numerical 
strong-disorder renormalization group to the dissipative 
transverse-field Ising model. ^^'^^ Their computer simula- 
tion results supported the above smeared-transition sce- 
nario but focused on the infinite-randomness physics aris- 
ing at intermediate energies. 

In this paper, we derive a comprehensive analytical ap- 
proach to the quantum phase transition of the transverse- 



field Ising chain coupled to dissipative baths of harmonic 
oscillators. The method is a generalization of Fisher's 
solution^"'^^ of the dissipationless case which treats on 
equal footing the effects of dissipation and disorder. For 
Ohmic damping, the theory captures the full crossover 
between the infinite-randomness (and quantum Griffiths) 
physics at higher energies where the dissipation is less 
important and the dissipation-dominated smeared tran- 
sition at low energies. The effects of sub-Ohmic dissipa- 
tion are qualitatively similar and also result in a smeared 
transition. In contrast, super-Ohmic dissipation is irrel- 
evant, and the transition is identical to that of the dissi- 
pationless chain. A short account of part of these results 
was already published in Ref. 18. 

Our paper is organized as follows: We define the 
model in Sec. II. In Sees. Ill and IV, we derive the 
renormalization-group recursion relations and the corre- 
sponding flow equations. They are solved in Sec. V. Sec- 
tion VI is devoted to phase diagrams, observables, and 
crossover effects. In Sec. VII, we discuss the applicabil- 
ity of our method for weak disorder, generalizations to 
higher dimensions, and the relevance of our results for 
experiments. We conclude and compare with related re- 
sults in the literature in Sec. VIII. 



II. THE HAMILTONIAN 

Our model Hamiltonian consists of a one-dimensional 
random transverse-field Ising model coupled to a bosonic 
bath/*^ 



H = Hi+He 



He 



Here, 



Hr 



Z2 ■^^^i'^i+l ^ zl ^"-^i 



(1) 



(2) 



is the Hamiltonian of the (dissipationless) random 
transverse-field Ising chain. The spin-1/2 degree of free- 
dom at site i is described by the Pauli matrices cr^, the 
local ferromagnetic interaction is represented by Ji, and 
the local transverse field hi controls the strength of the 
quantum fluctuations. 

Dissipation, i.e., damping of the quantum fluctuations, 
is introduced by coupling the z-component of each spin to 
an independent bath of quantum harmonic oscillators.^^ 
The bath Hamiltonian reads as 



Hf 



k.i 



t^fe/ 



at ,a,. 



(3) 



with ojfe^i being the frequency of the k-th harmonic oscil- 
lator coupled to the spin at site i, and a^. ■ (a J. ■) being 
the usual annihilation (creation) operators. (Note that 
we set h = 1.) Finally, the coupling between the spins 
and their respective baths is given by 



Hc^Yl ""^ H ^^^^ (" 



t 

kA 



"kA 



(4) 



with Afe^i being the strength of the interaction. 

All relevant information about the bosonic baths is en- 
coded in their spectral densities £i, which can be param- 
eterized as 

Ei{uj) = TT^ A^, ,,(5 (w - ujua) = -a^Lul~''uj'' {lj < Wc.j). 

k 

(5) 
Here, ujc^i is an ultraviolet cutoff energy, and ai is a di- 
mensionless measure of the dissipation strength. The mi- 
croscopic energy scale ujoa is often identifled with the 
cutoff iL)c,i- However, we need to keep the two ener- 
gies separate because the cutoff ujc,i will flow within our 
renormalization-group procedure, while ujQi (the bare 
cutoff energy) will not. Depending on the value of the 
exponent s, different types of dissipation can be distin- 
guished, s = 1 corresponds to the experimentally impor- 
tant Ohmic case. For s > 1, the bath is super-Ohmic, 
which leads to qualitatively weaker dissipation effects as 
there are fewer bath states at low energies. For s < 1, 
the bath is dubbed sub-Ohmic. 

We emphasize that in the model (1), each spin is 
coupled to its own (local) oscillator bath. Therefore, 
these baths damp the local order-parameter fluctuations, 
i.e.; they induce a long-range interaction in time. How- 
ever, they do not produce interactions between different 
sites. If we coupled all spins to the same (global) os- 
cillator bath, we would also obtain an effective interac- 
tion, mediated by the bath oscillators, between the spins 
at different sites. This interaction would be analogous 
to the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action between localized magnetic moments in metals. 

Quenched disorder is implemented by taking all pa- 
rameters in the Hamiltonian (1), i.e., the interactions Ji, 
the transverse fields hi, the dissipation strengths a^, as 
well as the bath energy scales uJca a-nd wq.^ to be in- 
dependent random variables. Actually, making one of 
these quantities random is sufficient because, under the 
renormalization-group flow (discussed in Sees. HI and 
IV) , the other quantities acquire randomness even if their 
bare values are not random. In the following, we thus as- 
sume that in the bare system only Ji and hi are random, 
with probability distributions Pi{J) and Rjih), while 
Q!i = a/ = a, uoA = ^ii and oJca = Wc are uniform. 
We also restrict ourselves to the experimentally most in- 
teresting case^^ of the bath energy w/ being the largest 
energy scale in the problem, loi S> hi^Ji. 



III. RENORMALIZATION GROUP 
RECURSION RELATIONS 

We now turn to the derivation of our theory. Our 
intent is to apply a real-space-based strong-disorder 
renormalization-group method to the Hamiltonian (1). 
This technique was introduced to tackle random antifer- 
romagnetic spin-1/2 chains'^"-^^ and has been successfully 
generalized to many disordered systems (for a review, see 
Ref. 22). Its philosophy is to successively integrate out 



the degrees of freedom with the highest local energies. 
This approach is justified a posteriori if the renormal- 
ized disorder strength becomes very large. 



A. Review of the dissipationless case 

The strong-disorder renormalization group of the dis- 
sipationless random transverse-field Ising chain (2) was 
derived and solved by Fisher. ^"'^^ One starts by identi- 
fying the highest local energy scale, Q = max{/ii, J^}, in 
the system. 

If the largest energy is an interaction, say, J2 , the low- 
energy states are those in which erf and cr| are parallel: 

|-|--|-) and I ). Therefore, the cluster of sites 2 and 3 can 

be recast as a single effective spin-1/2 degree of freedom, 
the effective magnetic moment of which is /I = /Lt2 + Ats . 
The transverse fields ft,2 and /13 will lift the degeneracy 

between |-|--|-) and | ) in second order of perturbation 

theory, yielding an effective transverse field h = /i2^3/^2 
for the new spin. 

In contrast, when the largest energy is a transverse 
field, say /12, the spin at site 2 is rapidly fluctuating 
between the |-|-) and |— ) eigenstates (in z direction). 
It hence does not contribute to magnetization and can 
be decimated from the chain. The neighboring spins 
CTi and (73 become connected via an effective coupling 
J = Ji J2//12, also obtained in second order of perturba- 
tion theory. 

This summarizes the renormalization-group recursion 
relations for the dissipationless case. Note that the re- 
cursion relations are symmetric under exchanging h and 
J, reflecting the self-duality of the Hamiltonian (2). 



B. Generalization to the dissipative case 

The generalization of the strong-disorder renormaliza- 
tion group to the full Hamiltonian (I) is not unique be- 
cause the bosonic baths can be incorporated in different 
ways. 

One possibility is implemented in the numerical work 
of Refs. 16,17. Here, the renormalization-group step still 
consists of identifying the largest of all the interactions 
and transverse fields, fl = inax{ J^, hi}. Before decimat- 
ing the corresponding high-energy degree of freedom, one 
integrates out the high-energy bath oscillators of the in- 
volved sites. This modifies the transverse fields and thus 
changes the recursion relations. Furthermore, the bath 
cutoff uJc becomes site dependent. While this scheme ap- 
pears to work very well numerically, we found it hard 
to implement analytically as it yields nonanalytic flow 
equations for the distributions of J, h, and fi. 

This difficulty is overcome by integrating out the oscil- 
lator degrees of freedom together with the lattice modes 
(related to strong interactions and transverse fields). We 
thus define the renormalization-group energy scale as 



il ~ inax{ Ji, hi, U!c/p}, where p ^ I is a dimension- 
less constant, the importance of which will become clear 
below in the context of adiabatic renormalization. A 
renormalization-group step now consists of decimating 
either an interaction Ji = r2, a transverse field hi = Q, 
or an oscillator with frequency LJk,i = P^- Upon lowering 
the renormalization-group energy scale from fl to 17 — dfJ, 
the bath cutoff Wc is thus reduced by pdf2 at all sites. In 
this sense, the renormalization of the baths is global, and 
not local as in the scheme of Refs. 16,17. This allows us 
to treat cj,, as a non-random variable. 

As it turns out, this modified scheme is analytically 
tractable. In the following subsections, we derive the 
resulting recursion relations. 



C. Adiabatic renormalization of the bosonic baths 

In order to understand the effects of dissipation on 
the quantum fluctuations, we first consider a single spin 
coupled to a bosonic bath, i.e., the well-known spin-boson 
problem^^ given by 



-ffsl 



-/icr^- 



k 



Xk 



+^a;fc (ala,^ + 



(6) 

The idea of the adiabatic renormalization is to eliminate 
high-frequency oscillators in a Born-Oppenheimer type of 
approximation. Let us assume that all oscillators of fre- 
quencies greater than ph with p ^ 1 can instantaneously 
follow the spin. We now wish to integrate out these fast 
oscillators. To this end, consider the Hamiltonian 



uik>ph 



Ho = a' J2 ^fc (4 + «A,- j + 2^ ^k[ a^flfe + 2 



E 

u>k>ph 



,t. 



E -^ Uk + ^) , 



u)k>ph 



(7) 



with 6j, = a^, + {Xkl^k)<^~ ■ Here, the sums are over 
the fast oscillators only, and an unimportant constant 
has been dropped. Notice that the difference between Ofe 
and bk simply represents a shift of t he oscillat or's rest 
position of magnitude xok = — 0'^AfeA/2/(TOfcW^). (Recall 
that Xk\/2rrikUJk = a\ + Uj^, with ruk being the mass of 
the fc-th oscillator.) Using the translation operator, the 
shifted ground state of oscillator k can be represented as 



|0fc,,.)=e""-PM0O = cxp 



; Afc 



,t 



\0k) , 



(8) 

where |Ofc) denotes the ground state of the unperturbed 
oscillator. Finally, the ground state of Hq is doubly de- 
generate with eigenvectors 



l*H 




n 10^-^ 

0Jk>ph 



± 



' n 10^^^ 

uik>ph 



(9) 



The transverse-field term in tlie Hamiltonian, Hi ~ 
— ha^, lifts this degeneracy. The tunnel splitting, i.e., 
the tunneling rate of the spin a between |+) and |— ), 
is given by the energy difference E+ — E^ = 2/i' = 
(vI/+|iJi|vI/^)-(*_|iJi|vI/_) = 2/in.^>p„(0fe,+ |0fe,_). 
Therefore, we can replace the spin coupled to the fast 
oscillators by an effective free spin with a renormalized 
transverse field 



h' = h \[ (Ofc,+ |Ofc,_)=/icxp 
2 /■"<= £ (w) 

T^ Jph 



-2 E (^'/^ 

Uk>ph 



/lexp 



-Auj 



(10) 



which depends exponentially on the spectral density (5). 
Now. as long as h' remains below the cutoff energy uj'^ = 
ph of the remaining oscillators, this procedure can be 
iterated until convergence. 

The qualitative features of this adiabatic renormal- 
ization procedure depend on the character of the 
dissipation.^^ For super-Ohmic dissipation (s > 1). h' 
always converges to a finite value ft,* > 0, regardless the 
dissipation strength a. Thus, the spin effectively decou- 
ples from the bath and remains tunneling in the presence 
of dissipation. In contrast, h* vanishes for sub-Ohmic dis- 
sipation (0 < s < 1) implying that the damping localizes 
the spin (at least in the limit of small h) . In the Ohmic 
case, the behavior depends on the dissipation strength 
a. As long as a < 1, the transverse field converges to 
a nonzero value h* = /i(p/i/ajc)"^^^~"^ while it converges 
to h* = for a > 1. Based on the heuristic arguments 
of Ref. 15, we therefore expect a smeared quantum phase 
transition in the sub-Ohmic and Ohmic cases. 

Let us now use the adiabatic renormalization of the 
oscillators within the strong-disorder renormalization- 
group scheme of the full dissipative random transverse- 
field Ising chain (1), as outlined in Sec. IIIB. Our intent 
is to lower the renormalization-group energy scale from 
ri to fi — dfi. To do so, we need to integrate out all oscil- 
lators in the frequency range between p{U, — ATI) and pf2. 
Because the interaction terms in (1) are diagonal in the 
a^ basis, they can simply be incorporated in Hq. Thus, 
when integrating out the high-frequency oscillators in all 
baths, the transverse fields renormalize according to 



K = hi 



exp 



npil. 

-aiOJj~^ I dcj 10^^"^ 

J p^2~pdfl 



1 - ttj 



r2 \ " ^ dn 



n. 



n 



(11) 



Here, we have used the relation uji — pQi between the 
bare renormalization group energy scale fli and the bare 
oscillator cutoff frequency w/. Notice that the renormal- 
ized fields /i'j do not depend on the arbitrary constant 
p> 1. 

Equation (11) reveals a crucial difference between the 
super-Ohmic, Ohmic, and sub-Ohmic cases. In the 



super-Ohmic case, the term renormalizing the transverse 
field is suppressed by a factor {fl/iljY^^, which vanishes 
as the energy scale fl goes to zero. This suggests that 
super-Ohmic dissipation becomes unimportant at suffi- 
ciently low-energy scales. Our calculations in the follow- 
ing will show that this is indeed the case. In contrast, the 
term renormalizing the transverse field remains finite or 
even diverges with 57 — >■ in the Ohmic and sub-Ohmic 
cases, respectively. 



D. Decimating a site 

Let us now focus on integrating out a site, say site 2, 
because its transverse field ^2 is at the renormalization- 
group energy scale Q. It is important to realize that if 
/i2 = f7, the field actually has reached its fully converged 
value /12 with respect to the adiabatic renormalization 
discussed in the last subsection (because all remaining 
oscillators have frequencies below p/i2). Hence, the spin 
(72 has decoupled from its bath, and decimation of the 
site proceeds as in the undamped case.^^ We consider 
Hq = — /i2C2 £is the unperturbed Hamiltonian, while the 
interactions of (T2 with its neighbors are treated perturba- 
tively: Hi ~ —{Jicrf + J20'|)cr|. In second order of per- 
turbation theory, the low-lying spectrum of the cluster of 
sites 1, 2, and 3 can be represented by 7?eff = ^J<^i<^3 
with 



J1J2 



(12) 



This means that the strongly fiuctuating spin (72 is re- 
moved from the chain, while the neighboring spins are 
coupled via a renormalized weaker bond J. 



E. Decimating an interaction 

Finally, the last possible decimation happens when 
an interaction, say J2, is at the renormalization group 
energy scale fi. In this case, we consider the cluster of 
sites 2 and 3 as well as its coupling to the neighboring 
sites 1 and 4. The unperturbed part of the Hamiltonian 
reads as Hq = — J2fT|o'|, while Hi^ = — /i2tT2 — ^acrf 
and Hiz = — Jicrfcrl 



Jso-Jct 



El2<E, 



A. 



(«L 



3(73^^4 



well 



H 



Idii 



are treated as perturbations. As in the undamped 
case, the low-energy states are those in which erf and 
(t| are parallel. Therefore, the cluster of sites 2 and 3 
can be replaced by a single effective spin a, the states 
|-|-) and |— ) of which represent \a2(J3) — \++) and 

\o'2<^3) = I ), respectively. Treating the perturbations 

to lowest non- vanishing order (first order of perturbation 
theory for Hiz and Hidiss but second order for Hi^), 



yields the effective Hamiltonian 



(a) Q, — ujc/p 



k 



-fe l+^k[ aia^ + 2 ) • (^^) 



The renormahzed transverse field h is smaller than either 
of the original ones and given by 



h2h3 
J2 



(14) 



The magnetic moment associated with the renormahzed 
spin a is the sum of the moments of (7| and a^ , 



/i = ^2 +M3, 



(15) 



while the dissipative bath coupled to a has spectral den- 
sity 



S{uj) = -zcet^i '^w'' with a = a2 + a^. 



(16) 



This means the renormahzed dissipation strength of a 
cluster increases with its magnetic moment. We can thus 
parametrize the local dissipation strength a^ in terms of 
the magnetic moment fj,i as 



a^i 



(17) 



This has the advantage that we do not have to separately 
keep track of the probability distribution of the dissipa- 
tion strengths during the renormalization-group flow. 

We note in passing that one can also include -ffidiss ui 
the unperturbed Hamiltonian, as was done in Ref. 17. 
This leads to more complicated algebra but results in 
the same recursions (13)-(16) at the approximation level 
considered here. 




Figure 1: (Color online) Decimation scheme. The magnitude 
of the transverse fields, interactions and magnetic moments 
are depicted by the length of the arrows, the thickness of the 
lines, and the size of the dots, respectively. The bosonic baths 
are depicted by the array of springs and oscillators in which 
lighter (small) oscillators are faster. Finally, the magnitude 
of the coupling to the bath (or the local dissipation strength) 
is depicted by the thickness of the spring connecting the spin 
to the bath, (a) Global adiabatic renormalization of the fast 
oscillators, (b.i) decimation of a strong local transverse field, 
and (b.ii) decimation of a strong local interaction. 



moment are given by Eqs. (14) and (15). The spectral 
function of the bosonic bath coupling to the effective spin 
has the same exponent s as the original baths, but the 
dissipation strengths increases according to (16) [see Fig. 
l(b.ii)]. This last renormalization is taken into account 
automatically if one dehnes the local at as in Eq. (17). 



F. Summary of the decimation procedure 



IV. FLOW EQUATIONS 



The entire decimation procedure resulting from the 
steps outlined in Sees. IIIC to HIE can be summarized 
as follows (see Fig. 1 for a schematic illustration). 

In lowering the renormalization group energy scale 
from O to fi — dr2, we first integrate out all oscillators 
with frequencies in the range [pQ —pdfl.pfl] in adiabatic 
renormalization. This modifies all transverse fields ac- 
cording to (11) [see Fig. l(a)[. 

After this global decimation, we integrate out all 
transverse fields and interactions in the energy range 
[n — dil,n]. (As these decimations are local they do 
not interfere with each other.) In the case of decimating 
a transverse field, the corresponding site and bath are re- 
moved from the system, yielding an effective coupling be- 
tween the neighboring spins given by (12) [see Fig. l(b.i)[. 
When integrating out an interaction, the corresponding 
spins are replaced by a single effective one. The trans- 
verse field acting on this spin and its effective magnetic 



The strong-disorder renormalization group consists 
in iterating the decimation steps of Sec. III. and so 
reducing the renormalization-group energy scale il = 
max{ Ji, /li, ojc/p} to zero. In this process, the proba- 
bility distributions of the interactions, transverse fields, 
and magnetic moments change. In this section, we de- 
rive the renormalization-group flow equations for these 
distributions. In contrast to the dissipationless case, 
where the magnetic moment does not enter the recur- 
sions for the transverse fields and interactions, ^"'^^ our 
recursion relation (11) couples the magnetic moments fi 
and transverse fields h (via a^ — fita). In addition to 
the distribution P( J; ft) of the interactions, we therefore 
need to consider the joint distribution of transverse fields 
and magnetic moments R{h,iJL;Q,). The last argument of 
these distributions makes explicit the dependence on the 
renormalization-group energy scale Vt. Occasionally, we 
will also use the reduced distribution of the transverse 



fields, Rh{h; fl) = f^ R{h, fi; n)dn. 



of a site can be found as 



A. Flow of P{J- Q.) 

Tlie distribution of the interactions J only changes 
when a site (a strong transverse field) is decimated. Upon 
reducing the energy scale from il to fi — dJ7, P{J; fl) thus 
transforms almost as in the undamped case^^: 



p{j; n - dn) [1 - dn {p{n: n) + r',^, (fi; o))] = 

P{J; fi) - dfii?;,,(fi; fi) / dJidJ2P(Ji; fi)P(J2; fi) 

J1J2 



6{J-.h) + 5{J-J2)-6[J- 



fi 



(18) 



where R'f^, {h' ; il) is the distribution of transverse fields 
after integrating out the fast bath oscillators in the en- 
ergy interval [pfi — pdfi,pfi] according to Eq. (11). The 
term between the brackets on the left-hand side of (18) 
guarantees the normalization of P{J; fi— dfi) (decimating 
either a site or an interaction each reduces the number 
of interactions in the system by one) . The second term 
on the right-hand side implements the decimation of a 
site (transverse field), which happens with probability 
dfii?Jj, (fi; fi). The first and second delta functions in this 
term remove the neighboring couplings while the third 
one insert the renormalized one, as given by (12). The 
integral over Ji and J2, weighted by P(Ji; fi)P(J2; fi), 
sums over all possible values of neighboring couplings. 

The only difference to the undamped case stems from 
the fact that the decimation of a site is determined by 
the distribution i?^^, (/i';fi) of the transverse fields after 
integrating out the fast bath oscillators rather than the 
distribution of the unrenormalized fields Rh{h] fi). Using 
(11), the probability R'/^, (fi; fi)dfi of having a decimation 



(a) Before 



Mc 



dO. h 



M^ 



(b) After 



^GfrT^ 



dn h' 



Figure 2: (Color online) Schematic of the possible transverse- 
field decimations in the h vs. /j. plane, (a) Before the integra- 
tion of the fast oscillators, the hatched area represents all the 
potential fields to be decimated, (b) After integrating out the 
oscillators, some of these fields have been shifted below il — dQ 
and are not eligible to be decimated anymore, only those in 
the triangular shaded area still are. The dashed line h — i^ 
maps into the dotted line fj, = {Q — h')/{a{i^/fli)"~^dQ,) after 
the adiabatic-renormalization step (11). 



i?;,,(fi;fi)dfi= / d/i / dh'R'{h',n;^) 

Jo Jn-dn 

dn / d/ii?(/i, ^; fi) 

1 - ajL^^ (fi/fi/)""^j i?/,(fi; fi)dfi, (19) 



where /Iji = /o A'^(^j A*; ^)dM//o P(f^, m; ^)d/i is the 
mean magnetic moment to be decimated at the energy 
scale fi. Equation (19) can be understood as follows (see 
Fig. 2 for an illustration) . Before integrating out the fast 
oscillators, the potential fields to be decimated are those 
in the range [fi— dfi, fi] highlighted in Fig. 2(a). However, 
after eliminating these oscillators, some of those fields are 
renormalized downward out of decimating region. Only 
the fields in the hatched area of Fig. 2(b) may still be 
decimated. 

Equation (19) implies that any cluster with magnetic 
moment larger than /Xc(fi) = l/(a(fi/fi/)''~"'^) will never 
be decimated. For Ohmic dissipation, s = 1, this means 
that there is an energy-independent cutoff cluster size 
He = 1/a (and thus a finite length scale). All clusters 
with fi > fic will become frozen as their transverse fields 
h iterate to zero with fi — )- 0. This will lead to the 
smearing of the quantum phase transition. In contrast, 
for super-Ohmic dissipation, s > 1, the cutoff cluster size 
/ic (fi) diverges with fi — )■ 0, which reflects the fact that 
super-Ohmic dissipation cannot prevent the tunneling of 
even the largest clusters. 



B. F\ow of R{h,n;n) 



The joint distribution R(h, fi; fi) of transverse fields 
and magnetic moments changes when an interaction is 
decimated. It also changes in response to the adiabatic 
renormalization of the fast oscillators. As a result, the 
fiow equation of R{h, fi; fi) is more complicated than in 
the dissipationless case:^^ 



R {h, fi;n- dn) [1 - dn {P{n; O) + R',^, {Q: O))] = 

R {h, fi; n) - dflPiCt; il) I dh2dn2dhsdn3R{h2, fi2-n)R{h3, fis; il) 



S {h — /i2) 5 [pL — ^2) + 5 {h— /13) S (/.( — /is) — S I h — 



n 



S (p - /i2 - Ms) 



d/i2d/X2i?(/i2, Ai2; ^) [S {h - ft.2) S (p. - p.2) - 6 {h - h'^) 5 {p.- ^2)] , 



(20) 



where h'2 is given by Eq. (11). The second term on the right-hand side of (20) implements the change in the fields 
and magnetic moments when an interaction is decimated. The third term implements the adiabatic renormalization 
due the decimation of the fast bath oscillators. The normalization term [between the brackets on the left-hand side 
of (20)] takes the same form as in (18) because decimating either a site or an interaction each reduces the number of 
sites (i.e., transverse fields) in the system by one. 

After some tedious but straightforward algebra, Eqs. (18) and (20) simplify to 
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(21) 



l + h—]Rih,ti;n) 



P{n; n) I dh2dii2Ahsdii^R{h2,ii2\ n)R{h3, /ig; n)S fh 



n 



5{iJ.- H2 - Ps)- 



(22) 



All terms proportional to a are due to dissipation. Those 
multiplying /Iq correspond to a reduced probability of 
decimating a site because the transverse fields are renor- 
malized downward by the dissipative baths. The addi- 
tional term (the second one on the right-hand side) in 
Eq. (22) corresponds to the global change of all trans- 
verse fields when the oscillators are integrated out. The 
flow equations of the undamped case^*''^"'^ are recovered 
as expected if we set a = 0. 

Moreover, in the undamped case, the magnetic mo- 
ments in (22) can be integrated out, resulting in a flow 
equations for the reduced distribution Rh{h;il) of the 
transverse fields. This is no longer the case for a 7^ 0, 
and the fixed-point distributions for fields and magnetic 
moments have to be obtained jointly. 



A. Logarithmic variables 



We introduce a logarithmic measure of the 
renormalization-group energy scale. 



T = in{ni/n) 



as well as the logarithmic variables 



(23) 



C = In (n/J) and /? = In (n/h) , 



(24) 



V. FORMAL SOLUTION OF THE FLOW 
EQUATIONS 

In this section, we solve the renormalization-group flow 
equations (21) and (22). To this end, we first change to 
logarithmic energy variables. We then perform a general 
rescaling transformation which allows us to identify fixed- 
point distributions of the flow equations. 



describing the interactions and transverse fields, 
respectively.^'^'^^ This is advantageous because the mul- 
tiplicative recursion relations (12) and (14) turn into ad- 
ditive ones: (^ = ^^ -|- (^2 and /3 = /32 -l- /Sa . 

The probability distributions tt and p of the logarith- 
mic variables can be defined as tt{(; r)dC, ~ P{J', ^)dJ, 
p(/3, /x; r)d/3d/i — R{h, fi;D,)dhdfi. Rewriting the flow 
equations (21) and (22) in terms of these distributions 



TTo - ( 1 - afi^e ''^ ^''^ I pp^o 



HO 



yields 

dn Ott 

+ (^1 - aPoe-^''^^'^) p^^oTT 4 TT, (25) 

dp dp 



dV 0/3 



C 
1 - ajL^e^^''^'^^^] ppfl - TTo 

"/^e ^' '^' -^+t^qP ® P, 
op 



piP.lJ) 



(26) 



where JIq = Jip^Q ~ Mo is the mean magnetic m.oment to 
be decimated, ttq = 7r(0; F) determines the probabihty 
of decimating an interaction while p^^o = /p(OiM;r)d/i 
determines the probability for decimating a transverse 

field. The symbol (X) denotes the convolution, tt (E) tt = 

/dCidC27r(Ci)7r(C2)^(C-Ci-C2). 

The origin of each term in (25) and (26) can be easily 
tracked. The first term on the right-hand side of each 
equation is due to the change of the variable (^ or /? when 
r changes. The second terms ensure the normalizations 
of the distributions. The last ones are responsible for 
implementing the recursion relations (12) and (14), and 
the third term in (26) implements the damping of the 
transverse fields according (11). 



B. Rescaling 

We now look for fixed-point solutions of the flow equa- 
tions (25) and (26), i.e., for distributions n and p that are 
stationary under the renormalization-group flow. From 
the solution of the dissipationless case,^°'^^ it is known 
that such flxed points only emerge after the variables ^ 
and /3 are appropriately rescaled. Let us consider the 
general transformations rj = C/fci^)^ ^ — /5///3(r), and 
1/ = p/ff^(r), where the scale factors /c,/3,,i(F) are func- 
tions of the logarithmic energy cutoff F only. Trans- 
forming the flow equations to the distributions V{rj) and 
Ti-{9, v) of the rescaled variables yields 
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where / = d//dF, j/q 
/7^(0,l^;F)d^^. 



Fe=o, 7^0=^(0; F), and Te^^c 



C. Fixed-point distributions 

The low-energy solutions of the flow equations describe 
stable phases or critical points. Before analyzing the 
equations in detail, we can already identify two simple so- 
lutions that correspond to the conventional paramagnetic 
and ferromagnetic phases. The key feature is a separa- 
tion of scales between the interactions and the transverse 
flelds. 

When all interactions J are larger than all transverse 
flelds h (i.e., their distributions do not overlap), the 
renormalization-group procedure will only decimate in- 
teractions until every site is included in the infinite clus- 
ter formed. This happens at a finite energy scale, viz., 
the minimum value of J in the chain. The system is thus 
in a conventional, homogenously ordered, ferromagnetic 
phase. 

Analogously, if all interactions are smaller than all 
transverse fields, only sites will be decimated under the 
renormalization group. Strictly, the interactions need to 
be smaller than all fully renormalized transverse fields in 
the sense of the adiabatic renormalization of Sec. Ill C 
(which also requires sufficiently weak dissipation). This 
is the conventional paramagnetic phase. As above, there 
is a finite characteristic energy, viz., the minimum value 
of the (renormalized) transverse fields in the system. 

In the absence of dissipation, these two conventional 
solutions are gapped, with the energy gap given by the 
final renormalization-group energy scale discussed above. 
In the presence of the dissipative baths, a true energy gap 
does not exist, but the final renormalization-group scale 
sets a characteristic energy for the magnetic excitations 
in the system. 

After these simple solutions, we now turn our atten- 
tion to the nontrivial ones at which the interactions J and 
transverse fields h compete with each other. This means, 
the distributions of the interactions and (renormalized) 
transverse fields overlap. In these cases, we expect the 
characteristic renormalization-group energy scale to van- 
ish, either because the system is critical or due to a 
quantum-Griffiths mechanism. We therefore take F — )> oo 
and search for F-independent solutions of the rescaled 
flow equations (27) and (28). In this analysis, we need 
to distinguish the different types of dissipation. 



1. Review of the dissipationless case 

In the dissipationless case,^^'^^ the analysis is simpli- 
fled by integrating out the magnetic moments. Their 
distribution can be found separately after the flxed-point 
distributions P*{t]) and 'R*g{9) have been obtained. 

There are three types of fixed points, (i) A critical so- 
lution for which duality requires that fc, = fp- In order 
to have a physical (normalizable and non-oscillatory) dis- 
tribution V{r]), we conclude from (27) that f(/f( must 
scale like 1//,^ with F, implying /^ ~ F.^^ Further anal- 
ysis (requiring that Vq > 0) leads to /^ ^ F"^, with 



= (1 + •\/5)/2 being the golden mean. The fixed-point 
distributions of the rescaled interactions and transverse 
fields are V*{x) = 'R*g{x) = e^^. There does not seem 
to be a simple closed-form expression for the joint fixed- 
point distribution of fields and moments, TZ*{9, v). 

(ii) On the ordered side of the transition, there is 
a line of fixed points parameterized by Vq describing 
the ferromagnetic quantum Griffiths phase. The scale 
factors are /^ = 1 and /^ = /^ = cxp(7'or). The 
fixed-point distributions are given by 'P*{ri) = Voe"^"^, 
n* (9, v) = 7^^,o7^e,oe-^''•°^l5 (7^e,o6l - 7^^,o^^), with 7^e,o 
and TZvfi being nonuniversal constants (see Appendix A) . 

(iii) Finally, the third type of solution is a line of fixed 
points on the disordered side of the transition. They 
are parameterized by TZq and describe the paramagnetic 
quantum Griffiths phase. The scale factors are fp — 1 
and Jq ~ exp(7?.or). The resulting fixed-point distri- 
butions for the interactions and transverse fields read 
V*{tj) = Voe-'^°'^ and 7^;(6') = 7^oe-^«'' with Vo a 
nonuniversal constant. The mean value of the magnetic 
moments scales as 7I ~ T,^^ which suggests that /^ = F. 
However, there is no fixed-point distribution R*{v)?^ In 
fact, if one assumes that TZ*{9,v) exists, one arrives at 
the incorrect conclusion that Ji = const. ^^ 



as the renormalization-group scale is increased from F to 
F-|-dF. As a result, f^/fp can diverge in the limit F — > 00, 
but never faster than 0(F). Moreover, in the presence of 
dissipation, integrating out the oscillators via Eq. (11) 
leads to a downward renormalization of the transverse 
fields and thus an increase of j3. A fixed-point solu- 
tion then requires a corresponding increase in fp. It is 
thus clear, that for any physical fixed-point solution of 
the flow equations (27) and (28), the dissipative terms, 
which are proportional to a.f^ji,e~^~^'^/fp, remain expo- 
nentially suppressed for any s > 1. 

How do we interpret these results physically? The cor- 
rections to the effective transverse fields due to super- 
Ohmic damping become smaller and smaller with de- 
creasing renormalization-group energy scale fJ and in- 
creasing size (magnetic moment) of the clusters. This 
holds both at criticality and in the surrounding quantum 
Griffiths phases. The quantum phase transition of the 
super-Ohmic random transverse-field Ising chain is thus 
identical to that of the dissipationless chain. Even du- 
ality should be recovered at criticality. All these results 
are in agreement with numerical work in Ref. 17. 



S. Ohmic dissipation 



2. Super-Ohmic dissipation 

To investigate the rescaled flow equations (27) and (28) 
in the super-Ohmic (s > 1) case, we first consider the 
stability of the dissipationless fixed points of Sec. V C 1 
with respect to the dissipative terms. In all dissipative 
terms in (27) and (28), the dissipation strength a appears 
in the combination af^e~^'^~^'^^/fi3. To analyze the be- 
havior of these terms close to the dissipationless fixed 
points, we use the scale factors /^ and fp discussed in 
Sec. VCl. For all three types of dissipationless fixed 
points (critical, ferromagnetic Griffiths, and paramag- 
netic Griffiths), we find that the dissipative terms are 
subleading and vanish in the limit F — > od for any s > 1. 
This means that the dissipationless fixed points are sta- 
ble against weak super-Ohmic damping. In other words, 
super-Ohmic dissipation is an irrelevant perturbation (in 
the renormalization-group sense) of the dissipationless 
random transverse-field Ising chain. 

We have also checked that the dissipationless fixed- 
point solutions are the only physical ones by tediously 
inspecting the flow equations (27) and (28). If we assume 
that the dissipative terms are the leading ones, we always 
arrive at unphysical fixed-point distributions. 

Instead of reproducing these calculations, we now 
present a more intuitive argument. /? and /i are both 
renormalized when an interaction is decimated, and both 
have the same additive recursion relations, /3 = /32 + /^a 
and jl = fj,2 + fJ'S- Therefore, one might expect that 
ffj, ^ fp- This is not entirely true because /3 depends 
explicitly on the value of the renormalization-group scale 
F and thus suffers an additional downward shift of — dF 



In contrast to the super-Ohmic case, Ohmic dissipation 
is not an irrelevant perturbation. The dissipative terms 
in the rescaled flow equations (27) and (28) are not expo- 
nentially suppressed and become important in all of the 
three regimes (critical, ferromagnetic Griffiths, paramag- 
netic Griffiths) . Hence, qualitatively different behavior is 
expected. 

In this section, we show that the Ohmic dissipation 
destroys both the paramagnetic quantum Griffiths phase 
and the quantum critical point. Only a ferromagnetic 
"Griffiths-type" phase survives. Before discussing the 
technical details, let us understand this result physi- 
cally. In a putative paramagnetic quantum Griffiths 
phase, there are arbitrarily large ferromagnetic clusters 
embedded in the paramagnetic bulk. As discussed at the 
end of Sec. IV A, clusters with magnetic moments larger 
than ^c = o;~^ can never be decimated in the case of 
Ohmic dissipation. In the low-energy limit, they there- 
fore become frozen and develop true static magnetic or- 
der. These static clusters will be aligned by any infinites- 
imally small interaction, giving rise to long-range order. 
The system is thus an inhomogeneous ferromagnet rather 
than a paramagnet. 

Let us now see how the same result follows from the 
fiow equations (27) and (28). We start by assuming 
that there is a paramagnetic quantum Griffiths phase. 
In such a phase, the majority of renormalization-group 
steps would be decimations of sites which renormalize 
the distribution V of the interactions. The convolu- 
tion term V ® V therefore has to be a leading term 
in the flow equation (27) to describe the downward 
flow of the interactions under repeated decimations of 
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sites. Moreover, as (1 — ca'Qf^)TZefl is the probability 
for decimating a transverse field, 1 — auofij. = c needs 
to remain nonzero and positive in the limit F — > oo. 
Thus, a physical solution from (27) requires that f^ di- 
verges exponentially and /^ is a constant which we can 
set to 1. Inserting this into (28) leads to the fixed- 
point solution TV = A{v) exp[— c7?,g o^/(l ~ ck'^/m)]' with 
A{i') being a function of v ensuring normalization and 
V* = 7'oCxp(-7'o^), with fp = 1, /^ = const, and 
Jq =exp(c7e^_or). 

However, this solution implies that there are no spin 
clusters with moments bigger than ^c — t^cffi = l/a. 
Thus, this solution cannot describe a quantum Griffiths 
phase in which, due to statistical disorder fluctuations, 
arbitrarily large magnetic clusters exist with small but 
nonzero probability. (Whenever the distributions of in- 
teractions and transverse fields overlap, it is always pos- 
sible to find an arbitrarily large region in which all inter- 
actions are bigger than all transverse fields.) 

Using the same arguments, it also follows that a quan- 
tum critical solution of the flow equations (27) and (28) 
(at which /^ has to diverge to reflect the diverging cor- 
relation length) can not exist. 

After having seen that Ohmic dissipation destroys both 
the paramagnetic quantum Griffiths phase and the quan- 
tum critical point (because no physical fixed-point so- 

lution is possible when the convolution term V ® V 
is present), we now consider ferromagnetic solutions. 
On the ferromagnetic side of the transition, almost all 
renormalization-group steps (at low energies) involve dec- 
imations of interactions, while almost no sites are deci- 

mated. Therefore, the convolution term V ®V drops out 
of the flow equation for F — >■ oo. Moreover, the flow of the 
interactions close to the fixed point should be analogous 
to that of the undamped case because their decimations 
are not influenced by the additional renormalization due 
to the bosonic baths. Therefore, we use the ansatz /^ = 1 
as in the undamped case. The flow equation (27) then 
yields 



V*{r,)=Voe 



'Vov 



(29) 



where Vo is a nonuniversal constant which parameterizes 
the line of ferromagnetic fixed points. 

Let us now turn our attention to the joint distribution 
TV {0, v) of transverse fields and magnetic moments at 
these fixed points. In the dissipationless case (Sec. V C 1) 
fields and moments rescale in the same way, /^ = /^g. In 
the presence of Ohmic dissipation, the transverse fields 
suffer additional downward renormalizations due to Eq. 
(11), which become more and more important as larger 
ferromagnetic clusters are formed with increasing F. We 
thus expect that f^/ fp — > with F — >■ oo. Under this as- 
sumption, the flow equation (28) greatly simplifles, and 
all dissipative terms drop out. This implies that the 
functional form of the fixed-point distribution in terms 
of 9 and v is identical to that of the undamped sys- 
tem in the ferromagnetic quantum Griffiths phase. As 



in Ref. 11. we can integrate over either v or j3 and ana- 
lyze the resulting reduced distributions TVg and 7?.* . This 
yields the fixed-point distributions TV^ = T^-eoc"^""^ and 
7?.* = Tl^oe~^"°'^ , just as in the undamped case. 

To fix the free parameters TZgQ and TZ^o 8-nd to find 
the joint distribution TZ*{6,i'), we need to take into ac- 
count the subleading term in the fiow equation (28) close 
to the fixed point. In contrast to the leading terms, the 
subleading ones do depend on the dissipation. The de- 
tails of this somewhat lengthy analysis are presented in 
Appendix A. We find /^ = e'^"^ and /^ = Fe'^"^ as well 
as TZ^o/TZgo = a. Moreover, the rescaled magnetic mo- 
ments and transverse fields become perfectly correlated 
at the fixed point, such that the joint distribution reads 
as 



n*il3,v) =noe-'^°''S{0-aiy) 



(30) 



The remaining nonuniversal constant TZq depends on the 
initial conditions. 

To summarize, instead of a quantum critical fixed point 
accompanied by lines of quantum Griffiths fixed points 
on each side of the transition, we only find a ferromag- 
netic solution that describes the physics for all overlap- 
ping distributions of interactions and renormalized trans- 
verse fields. To test the stability of this fixed point, 
we have computed-closed form solutions of the entire 
renormalization-group flow for some particularly simple 
initial distributions of transverse flelds and interactions 
(see Appendix A). In addition, we have implemented the 
recursion relations numerically^^ '^^ and verifled that the 
distributions flow toward our fixed point for many differ- 
ent initial distributions. 



4- Sub-Ohmic dissipation 

For the sub-Ohmic random transverse-field Ising chain, 
our strong-disorder renormalization group method can- 
not be applied because the adiabatic renormalization of 
the field does not capture all of the relevant physics of the 
problem. In the sub-Ohmic case, the adiabatic renormal- 
ization only works in the limit of weak transverse fields. ^^ 
In this limit, the dissipation suppresses the tunneling for 
any dissipation strength a. This means, the transverse 
fields all renormalize to zero for any a. Thus, within 
the adiabatic renormalization approach to the baths, the 
sub-Ohmic random transverse-field Ising chain does not 
have a quantum phase transition, and it is always in the 
ferromagnetic ground state. 

However, it is known from more advanced techniques 
that the sub-Ohmic spin-boson model with fixed weak 
dissipation does undergo a quantum phase transition at 
some nonzero value of the transverse field (see, e.g., Refs. 
27,28). In our context, this implies that small clusters 
fiuctuate but the dynamics of sufficiently large clusters 
freezes (their transverse fields vanish under renormaliza- 
tion) just like in the Ohmic case. We will briefiy discuss 
the resulting behavior in Sec. VII. 
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Figure 3: (Color online) Schematic phase diagrams in the 
transverse-field (/ityp) — temperature (T) plane, (a) Infinite- 
randomness critical-point scenario of the dissipationless and 
super-Ohmic chains. SO and SD denote the strongly or- 
dered and disordered conventional phases while WO and 
WD denote the ordered and disordered quantum Griffiths 
phases, (b) Smeared transition scenario for Ohmic dissipation 
(0 < Of <C 1), with the inhomogeneously ordered (lO) ferro- 
magnet replacing the WD Griffiths phase. The open circle 
marks the end of the tail of the smeared transition. The thin 
dashed lines represent finite-temperature crossovers between 
a quantum critical (QC) regime and a quantum paramagnet 
(PM) or ferromagnet (FM). The crossover temperatures T* 
and T** are discussed in the text. 



VI. PHASE DIAGRAM AND OBSERVABLES 

Using the fixed-point distributions found in Sec. V, 
we can now determine the phase diagram and compute 
thermodynamic observables close to the quantum phase 
transition. 



A. Phase diagram and crossovers 

Let us contrast the infinite-randomness critical-point 
scenario of the dissipationless random transverse-field 
Ising chain^*''^^ (which, according to Sec. VC2, also ap- 
plies to super-Ohmic damping) with the smeared transi- 
tion scenario emerging for Ohmic dissipation. 

The phase diagram of the dissipationless case as func- 
tion of temperature T and typical transverse field /ityp 
is sketched in Fig. 3(a) (keeping the typical interaction 
strength fixed). 

At zero temperature, the ferromagnetic-paramagnetic 



quantum phase transition is governed by a universal 
infinite-randomness fixed point (IRFP). It is accompa- 
nied on both sides by the gapless paramagnetic (weakly 
disordered) and ferromagnetic (weakly ordered) quan- 
tum Griffiths phases. For higher fields (to the right of 
the paramagnetic Griffiths phase), there is a conven- 
tional (strongly disordered) gapped quantum paramag- 
net where quenched disorder effects are irrelevant. For 
sufficiently weak fields (to the left of the ferromagnetic 
Griffiths phase), the system is a conventional (strongly 
ordered) gapped ferromagnet where quenched disorder is 
also irrelevant. 

As ferromagnetic order in one dimension is destabilized 
at any nonzero temperature, only crossovers will occur 
for T > [as depicted by the dashed lines in Fig. 3(a)]. 
Above tire critical point, there is tire quantum critical 
region which is characterized by activated scaling, i.e., 
the dynamical exponent is formally infinity. If the sys- 
tem is close to but not at the quantum critical point, it 
undergoes a crossover to one of the quantum Griffiths 
phases upon lowering the temperature. In these phases 
the dynamical scaling is conventional, i.e., the dynami- 



cal exponent is finite. 



l/|r|, where r 



Hyp 



the distance from criticality and he is the critical field. 
In the undamped case, duality requires the critical point 
to occur when the typical interactions equal the typical 
transverse fields, thus he = Jtyp- For super-Ohmic dis- 
sipation, the fields are renormalized downward by the 
oscillators. Using (10) we find he ~ Jtypexp{a/(s — 1)}. 
The crossover temperature T* is obtained from the acti- 
vated dynamical scaling at the IRFP, 



T* ~ exp(-const/|r|-'''^), 



(31) 



where :^ = 2 and ?/; = 1/2 are the correlation length ex- 
ponent and tunneling exponent of the IRFP, respectively, 
and the constant is of order unity. 

Having reviewed the undamped case, let us now discuss 
the case of Ohmic dissipation. For weak damping, a ^ 1, 
the renormalization group flow at high energies (where all 
clusters are still small) will be almost identical to that 
of the undamped system. Thus, the high temperature 
part of the phase diagram shown in Fig. 3(b) is analo- 
gous to its undamped counterpart. However, upon lower- 
ing the energy scale (i.e., the temperature), the clusters 
grow. In the undamped quantum Griffiths paramagnet 
discussed above their average magnetic moment grows as 
j[ ^ r'^'A(i-0) ln{ni/T), with = (l-f V5)/2, and at the 
quantum critical point it increases as /I ~^ [In(i7//T)]'^.^^ 
Therefore, it is clear that at some crossover tempera- 
ture T** , the cluster moment reaches the critical value 
fie — l/o^; E^nd dissipation becomes important. The value 
of T** can be estimated from the condition JL = 1/a; this 
gives 



rj^^ 



exp[ 



-const X {l/af/'*'] 



(32) 



very close the quantum critical point of the undamped 
system, and 
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rri^ 



cxp 



(-const X r^'^^'^-^V") 



(33) 



for larger transverse fields. 

At temperatures below T**, the quantum dynamics of 
the surviving clusters freezes completely (because their 
renormalized transverse fields vanish). These clusters 
thus behave like large classical spins. At zero temper- 
ature, they can be aligned by any infinitesimally weak 
interaction, leading to long-range order. Thus, the quan- 
tum Griffiths paramagnet (and quantum critical point) 
are replaced by an inhomogeneously ordered (10) ferro- 
magnetic phase [marked by the thick dashed line in Fig. 
3(b). This phase is due to local formation of magnetic or- 
der rather than a collective effect, which implies that the 
paramagnetic-ferromagnetic quantum phase transition is 
smeared. It also means that the end point of the tail of 
the smeared transition [marked by the open circle in Fig. 
3(b)[ is not a critical point. 



B. Observables 

We now turn our attention to observables. As shown 
in Sec. VC2, super-Ohmic dissipation is an irrelevant 
perturbation at the infinite-randomness critical point of 
the dissipationless random transverse-field Ising chain. 
Therefore, the low-energy behavior of observables in the 
super-Ohmic case is identical to that in the dissipation- 
less case, which is discussed in detail in Ref. 11 (see also 
Ref. 22 for a review of additional results) . This observa- 
tion agrees with numerical results of Schehr and Rieger.^^ 

Here, we therefore only consider the case of Ohmic 
dissipation. Moreover, we focus on the inhomogeneously 
ordered (10) ferromagnetic phase (i.e., on the tail of the 
smeared transition), which is the novel feature of the 
Ohmic chain. As we have seen in the analysis of the flow 
equations in Sec. V, weak dissipation does not change the 
behavior of the other phases qualitatively; it only leads 
to quantitative corrections to the undamped physics. 

To characterize the extent to which the bare (initial) 
distributions Pj and Rj of the interactions and transverse 
fields overlap, we introduce the probability of an interac- 
tion J being larger than a renormalized transverse field. 



w 



dJPiiJ) / dh,sRi{Ks), 



(34) 



where hcs = /i(/i/51/)"/'^~"-' is a local field fully renor- 
malized by the baths according to (11). If w = 1, all in- 
teractions are larger than all fields, and the system is in 
the conventional (strongly ordered) ferromagnetic phase. 
Conversely, if w = 0, all interactions are weaker than all 
transverse fields, putting the system in the conventional 
(strongly disordered) paramagnetic phase. Rare regions 
exist for < w < 1 with the inhomogeneously ordered 
regime (the tail of the smeared transition) corresponding 
to It; < 1. 



To leading approximation, the zero-temperature spon- 
taneous magnetization m is simply the magnetic mo- 
ment of the last spin cluster remaining after the 
renormalization-group energy scale 51 has been iterated 
to zero. In the tail of the smeared transition, u; ^ 1, this 
cluster is made up of a collection of well separated rare re- 
gions of minimum moment 1/a on which the interactions 
are larger than the transverse fields. The probability for 
finding a rare regions of 1/a sites is simply w'^^"'. We 
thus conclude that the magnetization in the tail of the 
smeared transition, u; ^ 1, behaves as 

A/a 



w 



(35) 



Note that this result implies that the dependence of the 
magnetization on the typical transverse field strength 
htyp is nonuniversal; it depends on the details of the 
bare distributions Pj and Rj. This nonuniversality 
agrees with heuristic results on other smeared phase 
transitions. ^^'^^"-^^ 

For weak dissipation and close to the location of the 
undamped quantum critical point, the magnetization can 
also be estimated by the fraction of undecimated spins 
when the undamped renormalization-group fiow reaches 
the crossover energy scale T** (because below T** al- 
most no sites will be decimated). This can be readily 
accomplished using the results of Ref. 11. The fraction 
of undecimated spins at scale T** is given roughly by 
the product of the density of surviving clusters at this 
scale n{T**) and their moment 'p{T**). As the moment 
is simply 1/a (from the definition of T**), we obtain 

m - a^/i'f'^')-^ (36) 

at the undamped quantum critical point, and 

m - a-^r" exp(-const x r^+'^'l'^'^-^^ /a) (37) 

on its paramagnetic side. On the ferromagnetic side, the 
magnetization is nonzero even in the absence of damping; 
weak dissipation thus provides a small correction only. 
The result (37) holds sufficiently close to the undamped 
quantum critical point. As the paramagnetic phase is 
approached for larger r, the behavior crosses over to (35) 
(see also Fig. 1 of Ref. 18). 

The (order-parameter) magnetic susceptibility x can 
be computed by performing the renormalization group 
until the energy scale il reaches the temperature T. All 
spins decimated in this process are rapidly fluctuating 
between up and down (i.e., pointing in the a:-direction 
due to their strong transverse field) and contribute little 
to the susceptibility. All clusters remaining at scale T 
are treated as free and contribute a Curie term ii^/T to 
the susceptibility. 

For temperatures above T**, the renormalization 
group fiow coincides with that of the undamped 
case. The temperature dependence of the susceptibil- 
ity is therefore identical to that of the dissipationless 
transverse-field Ising chain. ^^ In particular, in the quan- 
tum critical region, x increases as 

X - n{T)Ti\T)/T ^ M^i/T)Y^-''^/T, (38) 
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with vanishing T. In the paramagnetic Griffiths region, 
it behaves as 



^ ^ r'^+2'"^(l-'^)[ln(f77/T)]2T-l+l/^ 



(39) 



with z « l/(2r) being the dynamical exponent. 

Below T** , the renormalization group flow deviates 
from the undamped one, and the system crosses over into 
the inhomogeneously ordered (10) ferromagnetic phase. 
Here, the susceptibility can be calculated from the fixed- 
point solution of Sec. VC3. Following the same steps 
as in the ferromagnetic Griffiths phase of the undamped 
system, we obtain 



X 



T 



-l-l/z 



(40) 



with z' = l/Po being the low-energy dynamical exponent 
in the ferromagnetic phase [see Eq. (29)]. 

It is desirable to relate the dynamical exponent z' of 
the inhomogeneous ferromagnet to the bare distributions 
of J and h, or to the distance r from the undamped 
quantum critical point. The qualitative behavior is eas- 
ily discussed. On the ferromagnetic side of the undamped 
transition (r < 0), z' approximately agrees with its un- 
damped value because dissipation is a subleading correc- 
tion, hence, z' k, z k, — l/(2r). In the inhomogeneously 
ordered phase (i.e., on the paramagnetic side of the un- 
damped transition, r > 0), T'o is determined by the distri- 
bution of the interactions between the surviving clusters 
at energies below T**. With increasing transverse fields 
(increasing r), these clusters become rarer implying that 
the distribution of their interactions becomes broader. 
Thus, Vo decreases with increasing transverse fields in 
the tail of the smeared transition. In the far tail, w <^ 1, 
the clusters surviving at energy scale T** are essentially 
independent which means that their distances follow a 
Poisson distribution of width w^^/". This translates into 
the exponential distribution (29) of the coupling variables 
77, with Vo ~ w^'°' . Thus, the dynamical exponent 



1/Po 



'1/q 



(41) 



diverges at the endpoint of the tail where w = 0. 

Close to the undamped quantum critical point, the be- 
havior of z' can also be determined from the undamped 
renormalization group flow. We approximate the flow in 
the Griffiths paramagnet by disregarding the effects of 
dissipation until we reach the energy scale of T**. Below 
this scale, we assume that no further field decimations are 
performed. The value of the parameter Vq can then be 
obtained from the distribution V of the undamped Grif- 
fiths paramagnet at the energy scale F** = ln(ri//T**). 
This yields 



/ 1 
z ^ — exp const x 
r 



„Vi'(0-i)+i 



la] 



(42) 



Thus, z' increases with increasing r, in agreement with 
the above heuristic arguments. Our results for z' im- 
ply a dramatic change of the dynamical exponent when 



crossing over from the quantum Griffiths paramagnet to 
the inhomogeneously ordered ferromagnet at T** . Such 
behavior was indeed observed by means of a numerical 
implementation of the renormalization-group rules. ^^'^^ 



VII. DISCUSSION 
A. Applicability to weakly disordered systems 

In this section, we discuss to what extent our results 
apply to systems with weak or moderate bare disorder. 
The strong-disorder renormalization-group recursion re- 
lations derived in Sec. Ill become asymptotically exact 
in the limit of infinite disorder, but they are only ap- 
proximations for finite disorder. For this reason, the 
strong-disorder renormalization group needs to be com- 
plemented by other methods to discern the fate of weakly 
disordered systems. 

To analyze the limit of weak disorder, we can start 
from the Harris criterion,"^ which states that a clean crit- 
ical point is stable against weak disorder if its correla- 
tion length exponent v fulfills the inequality dv > 2. The 
clean Ohmic transverse-field Ising chain features a quan- 
tum critical point with an exponent value v « 0.638.'^^ 
Disorder is thus a relevant perturbation. Kirkpatrick 
and Belitz'^^ and Narayanan et al?'^'^^ studied a quan- 
tum order-parameter field theory with Ohmic dissipation 
by means of a conventional perturbative renormalization 
group. This method, which is controlled at weak dis- 
order, did not produce a stable critical fixed point, but 
runaway flow toward large disorder (where the strong- 
disorder renormalization group becomes better and bet- 
ter). This strongly suggests that our results apply for 
any amount of (bare) disorder. 

We now turn to the question of whether or not our 
results become asymptotically exact in the low-energy 
limit. To address this question, we need to discuss the 
widths of the distributions tt and p of the logarithmic 
interaction and transverse-field variables in the inhomo- 
geneously ordered ferromagnet. Using the fixed point 
distributions (29) and (30) and the corresponding scale 
factors /^ = 1 and fp ~ Fe''°'", we find the widths to be 
given by I/ttq = l/Vo and l/p^,o = Te^°^ /Uefi. Thus, 
the relative width of the transverse-field distribution di- 
verges, while the width of the interaction distribution 
remains finite in the limit F — >■ 00 (11 — >■ 0). 

This result can be easily understood by following the 
renormalization-group flow. Above the crossover energy 
scale Vt** (where the typical cluster size reaches 1/a), the 
fiow approximately coincides with that of a paramagnetic 
Griffiths phase of the undamped problem. Thus, (almost) 
only sites are decimated, which enormously broadens the 
distribution tt of the interactions. Below fi**, almost all 
decimations are interaction decimations. Therefore, tt 
does not renormalize further, and its width remains large 
but finite. In contrast, the transverse fields are driven to 
zero and their distribution broadens without limit. Our 
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results for the inhomogeneously ordered ferromagnet are 
thus not asymptotically exact, but represent a good ap- 
proximation. However, in the tail of the smeared tran- 
sition, w — ^ 0, the width of the interaction distribution 
diverges [see Eqs. (41) and (42)]. 

In this large-disorder limit, higher-order corrections 
to the recursion relations derived in Sec. Ill are sup- 
pressed. This includes couplings between different dissi- 
pative baths that appear upon integrating out a site. As 
these couplings only appear in fourth order of perturba- 
tion theory, they are very weak (and of short-range type) 
in the large-disorder limit. Thus, all corrections to our 
recursion relations are irrelevant in the renormalization- 
group sense, and our theory becomes asymptotically ex- 
act in the tail of the smeared transition. 

Finally, we discuss to what extent our results would 
change if the bare system had long-range interactions be- 
tween the different baths (or equivalently, a single global 
bath). In this case, the true zero-temperature behav- 
ior would be dominated by these interactions and differ 
from our results. However, sufficiently weak long-range 
couplings would become important only at very low tem- 
peratures [below even T**, see Fig. 3(b)]. Thus, our the- 
ory would remain valid in a broad temperature window. 
The precise value of the crossover energy scale depends 
on the bare value of the long-range interactions and on 
in the details of the model. It is a nonuniversal quantity 
and beyond the scope of this work. 



Sub-Ohmic dissipation 



As discussed in Sec. V C 4, our strong-disorder 
renormalization-group method can not be directly ap- 
plied to the sub-Ohmic case because there is no phase 
transition if the dissipative baths are treated within adi- 
abatic renormalization. However, a qualitative picture 
of the physics of the sub-Ohmic case can nonetheless be 
developed. 

The sub-Ohmic spin-boson model is known^^^^^ to un- 
dergo a continuous quantum phase transition from a 
fluctuating phase to a localized (frozen) phase as the 
transverse field is decreased (or the dissipation strength 
increased). Thus, the dynamics of sufficiently large 
spin clusters (which have small effective transverse fields 
and large dissipation strength) is always frozen. These 
frozen (classical) clusters can be aligned by any infinites- 
imal interaction, leading to magnetic long-range order. 
Just as in the Ohmic case, this mechanism replaces 
the paramagnetic Griffiths phase and the quantum crit- 
ical point by an inhomogeneously ordered ferromagnetic 
phase and smears the underlying ferromagnetic quantum 
phase transition. 



C. Higher dimensions 

We now discuss to what extent our results also ap- 
ply to dissipative random transverse-field Ising models 
in higher space dimensions. The renormalization-group 
steps of Sec. HI can be performed in complete analogy 
to the one-dimensional case, and the resulting recursion 
relations take the same form as Eqs. (11), (12), (14), 
and (15). However, decimating a site now generates ef- 
fective interactions between all pairs of its neighboring 
sites, while decimating an interaction leads to an effective 
site (cluster) that couples to a larger number of neigh- 
bors than before. Thus, both decimation steps change 
the connectivity of the lattice. Moreover, they introduce 
correlations between the remaining couplings. For this 
reason, an analytical treatment of the renormalization 
group flow appears to be impossible, even in the dissi- 
pationless case. However, the strong-disorder renormal- 
ization group has been implemented numerically in two 
and higher dimensions by keeping track of the irregu- 
lar "maze" of new interactions created by the decimation 
steps. ■^^'■^^ An analogous approach would be possible in 
the dissipative case. 

Even in the absence of a complete solution, the effects 
of dissipation on higher-dimensional random transverse- 
field Ising models can be understood qualitatively. The 
extra downward renormalization (11) of the transverse 
fields due to the baths is purely local, and thus also 
applies to higher dimensions. Just as in one dimen- 
sion, it reduces the probability of decimating a site, 
(1 — aJI^'i)Rh{^). Once the moment /x of a cluster 
reaches 1/a, it will never be decimated. We conclude 
that clusters with moments ^ > 1/a will freeze imply- 
ing that a critical fixed-point solution is impossible. As 
in one dimension, the quantum phase transition is there- 
fore smeared, and the paramagnetic Griffiths phase is 
replaced by an inhomogeneously ordered ferromagnetic 
phase. 

Although the zero-temperature quantum phase tran- 
sition of the higher-dimensional random transverse-field 
Ising models is similar to the one-dimensional case, their 
behavior at nonzero temperatures differs. In one dimen- 
sion, long-range order only exists at zero temperature 
(see our phase diagrams. Fig. 3). In higher dimensions, 
a long-range ordered ferromagnetic phase can exist at 
nonzero temperatures. It is important to note that the 
classical phase transition separating this phase from the 
paramagnet at nonzero temperatures is not smeared. At 
nonzero temperatures, even the largest clusters (rare re- 
gions) will thermally ffuctuate. Thus, a finite interaction 
is required to align them. This restores a sharp phase 
transition. ^^ 

In two or more space dimensions, disorder can also be 
introduced into the (dissipationless) transverse-field Ising 
model simply by diluting the lattice. For weak transverse 
fields, the ferromagnetic phase survives up to the percola- 
tion threshold p^ of the lattice. The quantum percolation 
phase transition emerging at pc shares many character- 
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istics with the infinite-randomness critical points of the 
generic random transverse-field Ising models.'^* Adding 
Ohmic dissipation suppresses the quantum fluctuations 
of sufficiently large percolation clusters. The results is an 
unusual classical super-paramagnetic phase. '^^ However, 
the percolation transition remains sharp as it is driven by 
the critical lattice geometry at the percolation threshold. 



D. Experiments 

To the best of our knowledge, an experimental real- 
ization of the microscopic Hamiltonian defined in Eqs. 
(1) to (4) has not been found yet. However, its order- 
parameter field theory, a one-component (j> theory with 
an inverse Gaussian propagator of the form G(q, w)"^ = 
r + q^ + jwl", also describes a number of experimentally 
important quantum phase transitions. Based on univer- 
sality, we expect the qualitative properties of these tran- 
sitions to be analogous to those of our model (1). 

For example. Hertz' theory'**' of the antiferromagnetic 
quantum phase transition of itinerant electrons leads to 
such an order-parameter field theory. In this case, the 
dissipation is Ohmic (s = 1); it is caused by the elec- 
tronic particle-hole excitations rather than abstract heat 
baths. For disordered itinerant ferromagnets, the dy- 
namic part of the Gaussian propagator takes the form 
|cj|/q^ rather than |a;|. Experiments on these disor- 
dered metallic quantum magnets often show unusual 
behavoir,'**"'*'^ and strong-disorder effects have been sug- 
gested as possible reasons. ''''* However, explicit verifica- 
tions of the quantum Griffiths and smeared-transition 
scenarios were missing for a long time. Only recently, 
some promising results have been found in several itiner- 
ant ferromagnets.'*'^"*^ For a more thorough discussion, 
the reader is referred to Refs. 48,49. 

We emphasize, that even though the leading terms in 
the order-parameter field theory of the dissipative ran- 
dom transverse- field Ising model (1) agree with the cor- 
responding terms in itinerant quantum magnets, their be- 
haviors are not completely identical. Importantly, in the 
Hamiltonian (1), each spin couples to its own indepen- 
dent dissipative bath. In contrast, in metallic magnets, 
all spins couple to the same Fermi sea. This produces 
an additional long-range interaction between the spins 
which has no counterpart in the Hamiltonian (1). It is 
mediated by the fermionic particle-hole excitations and 
known as the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction. 



tion. The quantum phase transition is thus in the same 
universality class as the dissipationless chain, i.e., it is 
governed by an infinite-randomness fixed point and ac- 
companied by quantum Griffiths singularities. 

In contrast, for Ohmic dissipation, the sharp quantum 
phase transition is destroyed by smearing because suffi- 
ciently large rare regions completely freeze at zero tem- 
perature. Note that this as a consequence of the interplay 
between disorder and dissipation as neither disorder nor 
dissipation alone can smear the quantum phase transi- 
tion of the transverse-field Ising chain. The behavior of 
the sub-Ohmic case is qualitatively similar to the Ohmic 
one: it also leads to a smeared transition. 

The results of this paper apply to the case of discrete 
Ising order-parameter symmetry. Systems with continu- 
ous 0{N) order-parameter symmetry behave differently. 
In contrast to sufficiently large Ising clusters, which freeze 
in the presence of Ohmic dissipation, 0{N) clusters of 
all sizes continue to fiuctuate with a rate that depends 
exponentially on their magnetic moment.^'' The result- 
ing sharp phase transition is controlled by an infinite- 
randomness critical point in the universality class of 
the undamped random transverse-field Ising model. ^^'^^ 
Super-Ohmic dissipation has weaker effects, ^^ and the 
renormalization-group flow is very similar to the dissipa- 
tionless case.^^ 

The difference between the discrete and continuous 
symmetry cases as well as the differences between the dif- 
ferent types of dissipation can all be understood in terms 
of a classification of rare-region effects according to the 
effective defect (rare- region) dimensionality**'^": If the 
rare-region dimensionality cIrr is below the lower criti- 
cal dimension d~ of the problem at hand, the transition 
is sharp and conventional. If df/R = d~ , the transition is 
still sharp but controlled by an infinite-randomness crit- 
ical point. Finally, if finite-size rare regions can order 
independently of each other {dim ^ d~), the transition 
is smeared. This classification applies to classical phase 
transitions and to quantum phase transitions that can 
be related to classical ones via the quantum-to-classical 
mapping. 

We expect our work or appropriate generalizations to 
shed light onto a variety of quantum phase transitions in 
disordered systems in which the order-parameter modes 
are coupled to additional noncritical soft modes. 



VIII. CONCLUSIONS 

To summarize, we have studied the quantum phase 
transition of the dissipative random transverse-field Ising 
chain via an analytical implementation of a strong- 
disorder renormalization-group method. We have shown 
that super-Ohmic dissipation is an irrelevant perturba- 
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Appendix A: The RG flow in the Ohmic 
ferromagnetic phase 

In this appendix, the flow equations (27) and (28) are 
studied in the ferromagnetic regime for Ohmic (s = 1) 
damping. 

For ferromagnetic fixed-point solutions, the convolu- 
tion term in (27) has to drop out, leading to 



fc ( dV\ 1 fdV 



(Al) 



As no sites are decimated, the typical value of C, cannot 
grow (as in the undamped case). Thus, /^ = 1 and 



V*{ii) = Voe 



-Von 



(A2) 



with Vq being an integration constant. With /^ = 1 and 
assuming that 



U/fp ^ 



(A3) 



(because we expect that fp grows faster than /^ due to 
the dissipation), the fixed-point equation (28) reads 



-f(-^f)-|(-^f 



Voin (E) n-n 



(A4) 



Interestingly, all dissipative terms (those involving a) 
drop out, and we recover the undamped fixed-point equa- 
tion. Therefore, we can integrate either over the fields 
or over the magnetic moments and analyze TZg or TZ^, 
separately. Physically relevant solutions require that all 
terms in (A4) survive in the limit F — > oo which implies 
fp/ fp = ffj./ ffi = const. Integrating over i' and Laplace 
transforming (A4), we arrive at 

= -cz^iV{z, 0) + 7'o7e*(z, 0) {il*{z, O) - l) , (A5) 

where il{v,n) = J^ e-''^-"'^7^(6', i/)d6'di/ and c > Vq is 
a constant. Equation (A5), which can be linearized via 
a* = — c((9iin;lnu)/7'o, then yields 



n*{v,0) 



1 



1 + (I'M) 



■Po/c 



(A6) 



where A is an integration constant. For c =/= Vq, the in- 
verse Laplace transformation gives TZg{0) ex 9^^^'^°/^ for 
6* > 1, which corresponds to Rl{h) ~ l/{h\\nh\^+'^"/'') 
for /i <C 1. Although this extremely singular distribu- 
tion is a fixed-point solution of the flow equations, it is 
not the attractive one for typical initial distributions. In- 
deed, this extreme fixed point can only be accessed if the 
hare distribution of fields already contains such strong 
singularity.^^ 



The less singular (and attractive) solutions correspond 
to c = Po- By inverse Laplace transforming Eq. (A6), 
we conclude that TZ^ ~ TZeoe~''^''°^ , which is identical 
to the undamped case.^^ The same analysis applies to 
TZl = TZuoe^^""" . Here, TZgo and TZ„o are nonuniver- 
sal constants. Moreover, since the fiow of the magnetic 
moments is identical to the undamped case, we conclude 



that /„ 



3'^°^ . Recall we cannot use the result of the 



undamped case in which ff; = /^ because of our assump- 
tion (A3). 

After finding the reduced fixed-point distributions TZg 
and TZ* , we now search for the joint fixed-point distribu- 
tion TZy, v) with jp/fp = fjf^ = Vo and f^/fp = 
in the limit F — > 00. The Laplace-transformed fiow equa- 
tion (A4) becomes 



o>7^ 07^ 

^ ^ ^"^'^ '^~^~ 

ov on 



TZ TZ-l 



(A7) 



which can be linearized via TZ = — 9„ilnu, where m = 
In w + In n. Hence, 



n*{v,n) = (l + /(w/n)V^) 



(A8) 



where /(fc) is a generic function constrained to diverge 
^ \nt/TZgo when fc — > 00, and to diverge ^ l/{\/kTZvi^) 
when fc — )- as dictated by the solution (A6). 

Our purpose now is to find the correct function /(fc) 
in (A6). It is easy to see that there are two limiting cases 
of the joint distribution TZ{v, n) that are compatible with 
the reduced distributions TZ*g and 7?.* found above, viz.. 



TZcoii{v, n) = [1 + — h 



V 

R0O 



n 
Rvo 



and 



'^uncorr i^i ^) 



-) 

ReoJ 



n 
Ruo 



(A9) 



(AlO) 



which, respectively, correspond to perfect correlations, 
TZcorA9,i^) = ReoRuoe-'^''°'^S{Rgoe - R^qi^), and no 
correlations, TZuncorr{0,v) = Reoe~^''°^Ruoe~^-°'' , be- 
tween transverse fields and magnetic moments. However, 
"T^uncorr docs not Satisfy the general form ( A8) . This sug- 
gests that correlations between the fields and magnetic 
moments arise from the fixed-point structure itself, re- 
gardless the strength of the damping a. Indeed, as we 
shall see below, the fiow is toward the perfect correlated 
case (A9). 

Since it is not possible to gather full information about 
f{v/n) by analyzing the fixed-point solution alone, we 
now analyze the renormalization-group fiow near the 
fixed point. This can only be accomplished if the F- 
dependence is not scaled out. We thus consider the un- 
sealed variables /3 and /z and use the fixed point (A8) as 



reference. Having in mind that /^ 



„Por 



, U/ff3 -> 0, 



and f[j/f[j — > Poj we keep only the main terms as well as 
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the leading corrections in the flow equation (26): 
dp dp , , 

9f ^^'^^^^ + ''^^P®P'Pl■ 



(All) 



It is easy to see that the neglected term corresponds to 
the third one on the right-hand side of (28), which is the 
one responsible for the normalization of 7?, due to a field 
decimation. The effects of damping thus enter only in 
the renormalization of the fields upon cooling the bath. 
The careful reader may notice that the neglected term 
is proportional to f^Jf/s as is the fourth one. Then, why 
is —aJlQPi3,op neglected and ~apdp/d/3 is not? The rea- 
son comes from the arising correlations between trans- 
verse fields and magnetic moments. Consider for instance 
the perfectly correlated case TZcon- With this choice. 



A*o = ft^v^ = /, 



/7?.corr(0, i^)dj/ 



= 0. 



This result is interesting because, as argued in Sec. IV, 
(1 — aJlf^pfjQjp is the probability of having a decimation 
of a field. Since the magnetic moments grow ~ /^ along 
the renormalization group flow, one might naively expect 
that the mean magnetic moment to be decimated, /Jqi 
also grows ~ /^. This would imply that the above prob- 
ability becomes negative (notice ppfi is nonzero). There- 
fore, correlations between the fields and moments ensure 
that this probability remains well-defined. 

By Laplace-transforming (All), we finally arrive at 



dp_ ^^fj9p 
dr dm 



7ro(p 



P)^ 



(A12) 



where p{b,m) = J e~^^ p' {13 , m)dl3 and p'{l3,m) = 
J e~'^^p{(3,p)dfi. (The term adp'{0,m)/dm also van- 
ishes because of the correlations between fields and mag- 
netic moments.) The first term on the right-hand side of 
(A12) is the desired leading correction to the fixed-point 
flow equation (A7). Inspection of (A12), guided by the 
fixed-point result (A8), yields 



p*{b,m;T) 



m \ ( m 



where F(y, b) is an integration constant which depends 
on b and on the combination y = T + ^, which implies 



m/b ex r, i.e., fp/f^ oc F. To see this clearly, let us 
compute F from some initial conditions. 

Consider firstly a perfectly correlated one, i.e.. 



-/I \ -> I b m 

p{b, m; 0) = Pcorr = H \ 

Ppo Pt^o 



(A13) 



In this case F{y,b) = 6(1 + appoy/p^o)e''°y /pj^o- Thus, 
p* (6, m; r) = p{b, to; F — )■ cxj) with 



p (6, to; F) = [l-\ ^ '^— 

as expected, and fp 



= e'"'^ 



(AM) 

Fe''«r(l + 



with ff, 

O(aF)^^). Notice that p* is of the perfectly correlated 
type (A9). Finally, we obtain the scale factor fp ~ F/^ 
which satisfies the requirements that f^/fp -^ and 
/m//m = ip/U = ^0 = TTo in the hmit F ^^ oo. 

To convince ourselves that p* is the relevant attractive 
fixed-point distribution, consider, for example, an initial 
condition which is totally at odds with (A13): the per- 
fectly uncorrelated case 



/5(6,to;0) 



In this case, F(j/, 6) = b[l 
which yields 



b 

P/30 



TO \ 



ay(P/30 + b)/pf,o]e^"y /pno 



p{b,m;T) = 1 



afpb + /pTO 



0{b,- 



In the limit F — > oo, and 5, to, — > 0, this distribution fiows 
to the attractive solution /5*(6,to,;F) (A14). Since both 
uncorrelated and perfectly correlated initial conditions 
flow to the same (correlated) fixed-point solution (A14), 
we expect all well-behaved (ferromagnetic) initial distri- 
butions to flow toward this fixed point. We have also 
checked this numerically in Ref. 26. After the inverse 
Laplace transformation, we finally obtain 

P*W,p;r) ^ ppoP^.oe-"^'>^6{ppoP-p^op), (A15) 

with p^o = Po/ffi and ppo = po/{afp). 
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